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We present measurements of Collins asymmetries in the inclusive process e"*"e“ —>■ hih 2 X, h\h 2 = 
KK, Ktt, tttt, at the center-of-mass energy of 10.6 GeV, using a data sample of 468 fb“^ collected by 
the BABAR experiment at the PEP-II B factory at SLAG National Accelerator Center. Considering 
hadrons in opposite thrust hemispheres of hadronic events, we observe clear azimuthal asymmetries 
in the ratio of unlike-sign to like-sign, and unlike-sign to all charged / 11/12 pairs, which increase with 
hadron energies. The Kn asymmetries are similar to those measured for the tttt pairs, whereas those 
measured for high-energy KK pairs are, in general, larger. 

PACS numbers: 13.66.Be, 13.87.Fh, 13.88.-|-e, 14.65.-q 


The Collins effect [H relates the transverse spin com¬ 
ponent of a fragmenting quark to the azimuthal distri¬ 
bution of final state hadrons about its flight direction. 
The chiral-odd, transverse momentum-dependent Collins 
fragmentation function (FF) provides a unique probe of 
quantum chromodynamics (QCD), such as factorization 
and evolution with the energy scale 

Additional interest has been sparked by the obser¬ 
vation of azimuthal asymmetries for pions and kaons 
in semi-inclusive deep inelastic scattering experiments 
(siDis) These are sensitive to the product of a 

Collins FF and a chiral-odd transversity parton distribu¬ 
tion function (PDF), one of the three fundamental PDFs 
needed to describe the spin content of the nucleon. Al¬ 
though these observations require nonzero Collins FFs, 
independent direct measurements of one of these chiral- 
odd functions are needed to determine each of them. 

In e“'"e“ annihilation, one can measure the product of 
two Collins FFs, and detailed measurements have been 
made for pairs of charged pions U No measure¬ 
ments are available for Kn and KK pairs, which are 
sensitive to different quark-flavor combinations, in par¬ 
ticular the contribution of the strange quark. Such mea¬ 
surements could be combined with SIDIS data to simulta¬ 
neously determine the Collins FFs and transversity PDF 
for up, down, and strange quarks 14 M- 

In this paper, we report the measurement of the Collins 
effect (or Collins asymmetry) for inclusive production of 
hadron pairs in the process e“'"e“ -T qq ^ hih2X, where 
hi, 2 = K^ or TT^, q stands for light quarks u or d or s, 
and X for any combination of additional hadrons. 

The probability that a transversely polarized quark 
(g^) with momentum direction k and spin fragments 


into a hadron h carrying zero intrinsic spin with momen¬ 
tum P/i, is defined in terms of unpolarized D\ and Collins 
fragmentation functions [ 2 ^ : 


(k X Vht) ■ S, 

zMh 


Df (z,P,t) = D\[z,PlT) + Ht\z.PlT) 

ZlVlh 

( 1 ) 

where Mh, PhT, and z = 2Ehl\fs are the hadron mass, 
momentum transverse to k, and fractional energy, re¬ 
spectively, with Eh its total energy and the e+e“ 
center-of-mass (c.m.) energy. The term including 11^ in¬ 
troduces an azimuthal modulation around the direction 
of the fragmenting quark, called Collins asymmetry. 

qq events, the q and g must be produced 


In e+e 


back-to-back in the e“'"e“ c.m. frame with their spin 
aligned. For unpolarized e"*" and e“ beams at BaBAR 
energies, the g and g spins are polarized along either the 
e"*" or e~ beam direction, so there is a large transverse 
component when the angle between the e+e“ and the qq 
axis is large. The direction is unknown for a given event, 
but the correlation can be exploited. Experimentally, the 
g and g directions are difficult to measure, but the event 
thrust axis n 21, 2^ approximates at leading order the 
qq axis, so an azimuthal correlation between two hadrons 
in opposite thrust hemispheres reflects the product of the 
two Collins functions. 

Figure [ 1 ] shows the thrust reference frame (RF12) [j^ . 
If not otherwise specified, all kinematic variables are de¬ 
fined in the e+e” c.m. frame. The Collins effect results in 
a cosine modulation of the azimuthal angle 4>i2 = 4 ‘iA 4>2 
of the di-hadron yields. Expressing the yield as a func¬ 
tion of 4>12 (after the integration over P^t), and dividing 
by the average bin content, we obtain the normalized 
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FIG. 1. (color online). Thrust reference frame (RF12). The 
azimuthal angles (f)! and 02 are the angles between the scatter¬ 
ing plane and the transverse hadron momenta Pti{t 2 ) around 
the thrust axis n. The polar angle 9th is the angle between h 
and the beam axis. Note that the difference between Pti(t 2 ) 
and PhT is that the latter is calculated with respect to the qq 
axis. 


rate 




.Ri2(</>12) = 1 + 


sin^ 6>th 
1 -I- cos^ 0th 


cos 012 • 




where the sum over the involved quark flavors is implied, 
0th is defined in Fig. [U zi( 2 ) is the fractional energy of the 
first (second) hadron, and the bar denotes the function 
for the q. Equation (I2|) involves only the moments of FF, 
which are defined as 



I^tI 

_ Mt . 


F{Zi,IkT^, 


( 3 ) 


with n = 0, 1, and Ikxl the transverse momentum of the 
quarks with respect to the hadrons they fragment into, 
which, in this frame, is related to the measurement of the 
transverse momenta of the two hadrons with respect to 
the thrust axis. 

Despite the simple form of the i?i 2 normalized rate, 
which involves only the product of moments of FFs, the 
RF12 frame comes with several downsides, among oth¬ 
ers of having to rely on Monte Carlo (MC) simulations 
when using the thrust axis as a proxy for the leading- 
order qq axis. An alternative frame is the analogue of the 
Gottfried-Jackson frame 23|, [2^ which uses the momen¬ 
tum of one hadron as a reference axis, and defines a single 
angle 0o between the plane containing the two hadron 
momenta and the plane defined by the beam and the ref¬ 
erence axis. We refer to this frame as RFO The 

corresponding normalized yield in the e+e“ c.m. system 


is 


23| 


in^ 0 

Ro(20o) =1 + 1 2 a cos200- 

1 + COS^ 02 

^ ^ _I (4) 

T[{2\x • kT h • PT - kT • PT)i^f Fi ] 

where 02 is the angle between the hadron used as ref¬ 
erence and the beam axis, h is the unit vector in the 
direction of the transverse momentum of the first hadron 
relative to the axis defined by the second hadron, and T 
is used to denote the convolution integral 


^ ^ f {pt + kx — qx) 

<? ( 5 ) 


with kx, Px, and qx the transverse momentum of the 
fragmenting quark, antiquark, and virtual photon from 
e“'"e“ annihilation, respectively, in the frame where the 
two hadrons are collinear, and X {X) = Di (Di) or 
Hi {Hi ). In this frame, specific assumptions on the 
kx-dependence of the involved functions are necessary 
to explicitly evaluate the convolution integrals. 


For this analysis we use a data sample of 468 fb~^ [25| 
collected at the c.m. energy ^/s ~ 10.6 GeV with the 
BaBAR detector [i^, at the SLAG National Accelera¬ 
tor Laboratory. We use tracks reconstructed in the silicon 
vertex detector and in the drift chamber (DGH) and iden¬ 
tified as pions or kaons in the DGH and in the Gherenkov 
ring imaging detector (DIRG). Detailed MG simulation is 
used to study detector effects and to estimate contribu¬ 
tion from various background sources. Hadronic events 
are generated using the Jetset [2^ package and undergo 
a full detector simulation based on Geant 4 (29l| . 

We make a tight selection of hadronic events in or¬ 
der to minimize biases due to detector acceptance and 
hard initial-state photon radiation (ISR), as they can in¬ 
troduce fake azimuthal modulations. Furthermore, final- 
state gluon (qqg) radiation also leads to angular asym¬ 
metries to be taken into account (^ . Requiring at least 
three charged tracks consistent with the e“'"e“ primary 
vertex and a total visible energy of the event in the lab¬ 
oratory frame Etot >11 GeV, we reject e+e“ —>• t+t“ 
and two-photon backgrounds, as well as ISR (qqg) events 
with the photon (one jet) along the beam line. About 
10% of ISR photons are within our detector acceptance, 
and we reject events with a photon candidate with energy 
above 2 GeV. We require an event thrust value T > 0.8 
to suppress qqg and BB events, and | cos 0th | < 0.6 so 
that most tracks are within the detector acceptance. 

We assign randomly the positive direction of the thrust 
axis, and divide each event into two hemispheres by 
the plane perpendicular to it. To ensure tracks are as¬ 
signed to the correct hemispheres, we require them to 
be within a 45° angle of the thrust axis and to have 
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z > 0.15. A “tight” identification algorithm is used to 
identify kaons (pions), which is about 80% (90%) efficient 
and has misidentification rates below 10% (5%). We se¬ 
lect those pions and kaons that lie within the DIRC ac¬ 
ceptance region with a polar angle in laboratory frame 
0.45 rad < 0iab < 2.46 rad. To minimize backgrounds, 
such as e+e” —?► followed by photon conversion, 

we require z < 0.9. 

We construct all the possible pairs of selected tracks 
reconstructed in opposite thrust hemispheres, and we cal¬ 
culate the corresponding azimuthal angles (pi, 4 ' 2 , and pn 
in the respective reference frames. In this way, we iden¬ 
tify three different samples of hadron pairs: KK, Ktt, 
and TTTT. To reduce low-energy gluon radiation and the 
contribution due to wrong hemispheres assignment, we 
require Qt < 3.5 GeV/c, where Qt is the transverse mo¬ 
mentum of the virtual photon from e“'"e“ annihilation in 
the frame where the two hadrons are collinear 23|. 

The analysis is performed in intervals of hadron frac¬ 
tional energies with the following boundaries: 0.15, 0.2, 
0.3, 0.5, 0.9, for a total of 16 two-dimensional (zi,Z 2 ) 
intervals. 

For each of the three samples, we evaluate the normal¬ 
ized yield distributions i?i 2 and Rq for unlike {U), like 
(L), and any charge combination (C) of hadron pairs as 
a function of (pi + (p 2 and 2(pQ, as shown in the left plot 
of Fig. [2] for KK pairs, for example. These combinations 
of charged hadrons contain different contributions of fa¬ 
vored and disfavored FFs, where a favored (disfavored) 
process refers to the production of a hadron for which one 
(none) of the valence quarks is of the same kind as the 
fragmenting quark. In particular, by selecting KK pairs, 
we are able to study the favored contribution of 

the strange quark, not accessible when considering tttt 
pairs only. 


1.04 
1.02 
b: 1 
0.98 
0.96 

- 3 - 2-10 1 2 3 - 3 - 2-10 1 2 3 

FIG. 2. (color online). Distributions of normalized yields 
(left plot) for unlike ([/), like (L), and any charge combination 
(C) of KK pairs, and their double ratios (right plot) in RF12. 



The normalized distributions can be parametrized with 
a cosine function: = ba + a), cos pic, where a = 0, 12 

indicates the reference frames, i = U, L, C the charge 
combination of hadron pairs, and /3i2(o) = <Pi 2 i 2 (po)- 
The distributions are strongly affected by instru¬ 


mental effects. In order to reduce the impact of the detec¬ 
tor acceptance, as well as any remaining effect from gluon 
bremsstrahlung [ 2 ^ , we construct two double ratios (DR) 
of normalized distributions, R^/Ra R^/R^. The 
two ratios give access to the same physical quantities as 
the independent Rp^, that is the favored and disfavored 
FFs, but in different combinations. We report the results 
for both kind of DRs, which are strongly correlated since 
they are obtained by using the same data set. These are 
shown in the right plot of Fig. [2] for KK pairs in RF12. 
At first order, the double ratios are still parametrized by 
a function that is linear in the cosine of the corresponding 
combination of azimuthal angles: 

i?j,^ = :§=3Bb+A^cos/3„, (6) 

Rol 

with B and A free parameters, and i, j = U, L, C. The 
constant term B must be consistent with unity, while 
A contains the information about the favored and disfa¬ 
vored Collins FFs. 

We fit the binned R]i distributions independently for 
KK, Ktt, and tttt hadron pairs. Using the MC sample, 
we evaluate the K/tt (mis)identification probabilities for 
the 16 (zi, Z 2 ) intervals in each of the three samples. For 
example, the probability that a true KK pair is re¬ 
constructed as KK pair is about 90% on average, slightly 
decreasing at higher momenta, while the probability 
that a true Ktt pair is identified as KK is about 10%, 
and is negligible. 

The presence of background processes could introduce 
azimuthal modulations not related to the Collins effect, 
and modifies the measured asymmetry as follows: 



with nm = KK, Ktt, tttt, and i = cc, BB, . In 

Eq. [3 Anm are the true Collins asymmetries produced 
from the fragmentation of light quarks in the three sam¬ 
ples, A\^ is the f-th background asymmetry contribu¬ 
tion, and are the fractions of reconstructed kaon 

pairs coming from uds and background events, calcu¬ 
lated from the respective MC samples. By construction, 
Fi + Fuds = 1. A similar expression holds for Ktt and 
TTTT samples. 

Previous studies [H show that e“''e —>■ BB and t+t 
events have negligible F^^ < 2%, and FV+t- sig¬ 
nificantly different from zero only for the tttt sample at 
high z values. Since Fee can be as large as 30%, and A‘^‘^ 
are unknown, we determine App^ in Eq. 0 from sam¬ 
ples enhanced in cc by requiring the reconstruction of 
at least one D*^ meson from the decay D*^ —>■ D^tt^, 
with the 1?° candidate reconstructed in the following 
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four Cabibbo-favored decay modes: K~7t~^, isr“7r+7r“7r+, 
Kg7r~^Tr~, and . These modes are assumed to 

provide a representative sample of tttt, Ktt, and KK 
pairs to be used in the correction, an assumption that 
is strengthened by the observation that the background 
asymmetries for those modes were found to be consis¬ 
tent. We solve the system of equations for A^^^, 

for the standard and charm-enhanced samples, 
and we extract simultaneously the Collins asymmetries 
Akk, Aktv, and corrected for the contributions of 
the background and K/tt (mis)identification. The domi¬ 
nant uncertainties related to this procedure come from 
the limited statistics of the I?*-enhanced sample and 
from the fractions Fi. The uncertainties on the frac¬ 
tions are evaluated by data-MC comparison and amount 
to a few percent. All these uncertainties are therefore 
included in the statistical error of the asymmetries ex¬ 
tracted from the system of Eq. ©• 

We test the DR method on the MC sample. Spin ef¬ 
fects are not simulated in MC, and so the DR distribu¬ 
tions should be uniform. However, when fitting the dis¬ 
tributions for reconstructed KK pairs with Eq. m , we 
measure a cosine term in the full sample of 0.004 ± 0.001 
and 0.007 ± 0.001 in the RE12 and REO frames, respec¬ 
tively, indicating a bias. Smaller values are obtained for 
Ktt and tttt pairs (see Fig. Studies performed on the 
MC samples, both at generation level and after full sim¬ 
ulation, demonstrate that the main source of this bias is 
due to the emission of ISR, which boosts the hadronic 
system and distorts the angular distribution of the final 
state particles, resulting in azimuthal modulations not re¬ 
lated to the Collins effect. This effect is more pronounced 
for KK pairs due to the lower multiplicity with respect 
to the other two combinations of hadrons. Assuming the 
bias, which is everywhere smaller than the asymmetries 
measured in the data sample in each bin, is additive, we 
subtract it from the background-corrected asymmetry. 

Using the uds MC sample, or light quark e“'"e“ —>■ qq 
MC events, we study the difference between measured 
and true azimuthal asymmetries. The asymmetry is in¬ 
troduced into the simulation by reweighting the events 
according to the distribution 1 ± a • cos^®®", where we 
use different values of a ranging from 0 to 8% with posi¬ 
tive (negative) sign for U {L and C ) hadron pairs, and 
01®" are the azimuthal angles combinations calculated 
with respect to the true qq axis in RF12, or the gen¬ 
erated hadron momentum in REO. The reconstructed 
asymmetries in RE 12 are systematically underestimated 
for the three samples of hadron pairs, as expected since 
we use the thrust axis instead of the qq axis, while they 
are consistent with the simulated ones in REO, where 
only particle identification and tracking reconstruction 
effects could introduce possible dilution. Since we mea¬ 
sure the same dilution for KK^ Ktt, and tttt samples, the 
asymmetry is corrected by rescaling Akk, Ak-w, and Atttt 
using the same correction factor, which ranges from 1.3 


to 2.3 increasing with z, as shown in Fig. [3l No correc¬ 
tions are needed for the asymmetries measured in REO. 
The uncertainties on the correction factors are assigned 
as systematic contributions. 



0.15 0.2 0.3 0.5 0.2 0.3 0.5 0 . 20.3 0.5 0.2 0.3 0.5 0.9 

0.9/0.15 0.9/0.15 0.9/0.15 


0.15 0.2 0.3 0.5 0.9 


FIG. 3. (color online). Correction factors for the dilution of 
the asymmetry due to the difference between the thrust and 
the qq axis. The open (full) markers, triangles and circles, 
refer to the U/L and U/C double ratios in the RF12 (RFO) 
frame, respectively. The 16 ( 21 , 22 ) bins are shown on the 
x-axis: in each interval between the dashed lines, 21 is chosen 
in the following ranges: [0.15,0.2], [0.2,0.3], [0.3,0.5], and 
[0.5, 0.9], while within each interval the points correspond to 
the four bins in 22 . 

All systematic effects, if not otherwise specified, are 
evaluated for each bin of 2 . The main contribution comes 
from the MC bias. We compare the bias results from 
the nominal selection, with those obtained by requiring 
different cuts on Atot, and/or by changing the detector 
acceptance region for the hadrons. The largest variation 
of the bias is combined in quadrature with the MC sta¬ 
tistical error and taken as systematic uncertainty. The 
effects due to the particle identification are evaluated us¬ 
ing tighter and looser selection criteria. The largest de¬ 
viations with respect to the nominal selection are taken 
as systematic uncertainties: the average relative uncer¬ 
tainties are around 10%, 7%, and 5% for the KK, Ktt, 
and TTTT pairs. Fitting the azimuthal distributions us¬ 
ing different bin sizes, we determine relative systematic 
uncertainties, which are not larger than 5%, 1.9%, and 
1% for the three samples. The systematic uncertainty 
due to the Etot cut is obtained by comparing the mea¬ 
sured asymmetries with those obtained with the looser 
selection Etot > 10 GeV. The average systematic con¬ 
tribution is around 10% for the three samples in both 
reference frames. We use different htting functions with 
additional higher harmonic terms. No significant changes 
in the value of the cosine moments with respect to the 
standard fits are found. As a cross-check of the double 
ratio method we ht the difference of K distributions, and 
we compare the two results. The difference between the 
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0.15 0.2 0.3 0.5 0.2 0.3 0.5 0 . 20.3 0.5 0.2 0.3 0.5 0.9 

0 . 9 / 0.15 0 . 9 / 0.15 0 . 9 / 0.15 


0.15 0.2 0.3 0.5 0.2 0.3 0.5 0 . 20.3 0.5 0.2 0.3 0.5 0.9 

0 . 9 / 0.15 0 . 9 / 0.15 0 . 9 / 0.15 


0.15 0.2 0.3 0.5 0.9 


0.15 0.2 0.3 0.5 0.9 


FIG. 4. (color online). Comparison of U/L (top) and U/C (bottom) Collins asymmetries in RF12 (left) and RFO (right) for 
KK, Kir, and tttt pairs. The statistical and systematic uncertainties are represented by the bars and the bands around the 
points, respectively. The 16 ( 2 : 1 , 2 : 2 ) bins are shown on the x-axis: in each interval between the dashed lines, 21 is chosen in the 
following ranges: [0.15, 0.2], [0.2, 0.3[, [0.3, 0.5[, and [0.5, 0.9], while within each interval the points correspond to the four bins 
in 22 . 


two procedures is negligible for Ktt and tttt pairs, while 
it reaches 1% and 3% for kaon pairs in RF12 and RFO, 
respectively. All the other systematic contributions are 
negligible 11 [. 

The Collins asymmetries measured for the 16 two- 
dimensional ( 21 , 22 ) bins, for reconstructed KK, Kn, and 
TTTT hadron pairs, are shown in Fig. 0] for RF12 and RFO, 
and are summarized in Tables mini and nil- The asym¬ 
metries are corrected for the background contributions 
and K/t: contamination following Eq. ([71), the MC bias 
is subtracted, and the corrections due to the dilution ef¬ 
fects are applied. The total systematic uncertainties are 
obtained by adding in quadrature the individual contri¬ 
butions, and are represented by the bands around the 
data points. 

An increasing asymmetry with increasing hadron ener¬ 
gies is visible for the U/L double ratio in both reference 
frames. The largest effects, but with less precision, are 
observed for KK pairs, for which AY 2 is consistent with 
zero at low 2 , and reaches 22% in the last 2 bin, while 
somewhat smaller values are seen for tttt and Kt: pairs. 
In particular, at low ( 21 , 22 ) bins for tttt pairs is 

nonzero, in agreement with the behavior observed in mi. 
The small differences between the two data sets are due 
to the different kinematic region selected after the cut on 
cos 9th- The A^’^ asymmetry is smaller than A^^ in all 
cases, and, for the KK pairs, the rise of the asymmetry 
with the hadron energies is not evident. We also note 
that the asymmetries for the KK pairs are larger than 
the others when the U/L ratio is considered, while they 
are at the same level, or lower, when they are extracted 


from the U/C ratio. 

In summary, we have studied for the first time in 
annihilation the Collins asymmetry for inclusive 


+0— 


e^e 


production of KK and Ktt pairs as a function of 
( 21 , 22 ) in two distinct reference frames. We measure 
the azimuthal modulation of the double ratios UjL and 
UjC, which are sensitive to the favored and disfavored 
Collins FFs for light quarks. We simultaneously extract 
also the Collins asymmetries for tttt pairs, which are 
found to be in agreement with those obtained in previous 
studies 11, 1^. The results reported in this paper and 
those obtained from SIDIS experiments can be used in 
a global analysis to extract the favored contribution of 
the strange quark, and to improve the knowledge on the 
u and d fragmentation processes dUlli. 
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KK sample 


Zl 

( 21 ) 

Z2 

{Z2) 

<sin^ Sth) 

tin 

-2 

) 

dVC tin 

,-2 

'l 

(1+cos^ e-fh.) 


^12 
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) 

[0.15,0.2] 

0.175 

[0.15,0.2] 

0.175 

0.797 

0.88 

± 

1.43 

± 

0.73 

0.62 

± 

1.02 

± 

0.51 

[0.15,0.2] 

0.175 

[0.2, 0.3] 

0.247 

0.794 

1.96 

± 

1.32 

± 

0.67 

0.88 

± 

0.92 

± 

0.43 

[0.15,0.2] 

0.175 

[0.3, 0.5] 

0.381 

0.794 

3.38 

± 

1.23 

± 

0.73 

1.08 

± 

0.81 

± 

0.43 

[0.15,0.2] 

0.175 

[0.5, 0.9] 

0.608 

0.786 

7.32 

± 

1.06 

± 

0.97 

2.06 

± 

0.72 

± 

0.59 

[0.2, 0.3] 

0.246 

[0.15,0.2] 

0.175 

0.794 

2.06 

± 

1.33 

± 

0.68 

0.95 

± 

0.93 

± 

0.44 

[0.2, 0.3] 

0.247 

[0.2, 0.3] 

0.247 

0.792 

3.78 

± 

1.22 

± 

0.58 

1.30 

± 

0.84 

± 

0.36 

[0.2, 0.3] 

0.247 

[0.3, 0.5] 

0.382 

0.792 

7.44 

± 

1.14 

± 

0.59 

2.06 

± 

0.72 

± 

0.35 

[0.2, 0.3] 

0.247 

[0.5, 0.9] 

0.608 

0.783 

10.91 

± 

0.98 

± 

0.91 

2.63 

± 

0.58 

± 

0.48 

[0.3, 0.5] 

0.381 

[0.15,0.2] 

0.175 

0.794 

5.34 

± 

1.28 

± 

0.74 

1.73 

± 

0.84 

± 

0.44 

[0.3, 0.5] 

0.381 

[0.2, 0.3] 

0.247 

0.792 

8.74 

± 

1.17 

± 

0.59 

2.45 

± 

0.74 

± 

0.35 

[0.3, 0.5] 

0.382 

[0.3, 0.5] 

0.382 

0.792 

10.97 

± 

1.31 

± 

0.63 

2.41 

± 

0.65 

± 

0.34 

[0.3, 0.5] 

0.383 

[0.5, 0.9] 

0.610 

0.784 

9.84 

± 

1.16 

± 

0.92 

1.71 

± 

0.52 

± 

0.45 

[0.5, 0.9] 

0.609 

[0.15,0.2] 

0.175 

0.785 

6.15 

± 

1.09 

± 

0.98 

1.78 

± 

0.74 

± 

0.60 

[0.5, 0.9] 

0.608 

[0.2, 0.3] 

0.248 

0.783 

11.75 

± 

1.03 

± 

0.91 

2.81 

± 

0.60 

± 

0.48 

[0.5, 0.9] 

0.610 

[0.3, 0.5] 

0.383 

0.784 

7.40 

± 

1.13 

± 

0.91 

1.11 

± 

0.52 

± 

0.45 

[0.5, 0.9] 

0.615 

[0.5, 0.9] 

0.615 

0.776 

22.36 

± 

2.09 

± 

1.69 

2.63 

± 

0.61 

± 

0.62 

Zl 

(zi) 

Z 2 

( 22 ) 

(sin^ 62 ) 

aVL /in 

-2 

'l 

aUC /in 

1-2 

) 

(l+cos^ 62 ) 


^0 



) 




[0.15,0.2] 

0.175 

[0.15,0.2] 

0.175 

0.739 

2.41 

± 

1.44 

± 

1.25 

0.82 

± 

1.02 

± 

0.58 

[0.15,0.2] 

0.175 

[0.2, 0.3] 

0.247 

0.736 

1.66 

± 

1.31 

± 

0.74 

0.53 

± 

0.92 

± 

0.40 

[0.15,0.2] 

0.175 

[0.3, 0.5] 

0.381 

0.750 

1.33 

± 

1.24 

± 

0.71 

0.23 

± 

0.81 

± 

0.39 

[0.15,0.2] 

0.175 

[0.5, 0.9] 

0.608 

0.751 

0.24 

± 

1.02 

± 

0.90 

-0.13 

± 

0.71 

± 

0.55 

[0.2, 0.3] 

0.246 

[0.15,0.2] 

0.175 

0.739 

1.95 

± 

1.32 

± 

0.75 

0.61 

± 

0.93 

± 

0.41 

[0.2, 0.3] 

0.247 

[0.2, 0.3] 

0.247 

0.736 

3.28 

± 

1.21 

± 

0.69 

1.00 

± 

0.84 

± 

0.33 

[0.2, 0.3] 

0.247 

[0.3, 0.5] 

0.382 

0.749 

3.69 

± 

1.14 

± 

0.67 

0.90 

± 

0.72 

± 

0.31 

[0.2, 0.3] 

0.247 

[0.5, 0.9] 

0.608 

0.750 

6.05 

± 

0.96 

± 

0.88 

1.49 

± 

0.58 

± 

0.44 

[0.3, 0.5] 

0.381 

[0.15,0.2] 

0.175 

0.738 

2.62 

± 

1.27 

± 

0.72 

0.67 

± 

0.84 

± 

0.39 

[0.3, 0.5] 

0.381 

[0.2, 0.3] 

0.247 

0.736 

4.76 

± 

1.17 

± 

0.67 

1.21 

± 

0.74 

± 

0.31 

[0.3, 0.5] 

0.382 

[0.3, 0.5] 

0.382 

0.749 

2.99 

± 

1.18 

± 

0.78 

0.73 

± 

0.63 

± 

0.31 

[0.3, 0.5] 

0.383 

[0.5, 0.9] 

0.610 

0.750 

4.18 

± 

1.06 

± 

0.92 

0.77 

± 

0.52 

± 

0.41 

[0.5, 0.9] 

0.609 

[0.15,0.2] 

0.175 

0.731 

3.30 

± 

1.10 

± 

0.91 

0.90 

± 

0.74 

± 

0.56 

[0.5, 0.9] 

0.608 

[0.2, 0.3] 

0.248 

0.728 

5.83 

± 

1.00 

± 

0.88 

1.42 

± 

0.60 

± 

0.44 

[0.5, 0.9] 

0.610 

[0.3, 0.5] 

0.383 

0.743 

2.52 

± 

1.04 

± 

0.92 

0.38 

± 

0.51 

± 

0.41 

[0.5, 0.9] 

0.615 

[0.5, 0.9] 

0.615 

0.743 

6.81 

± 

1.72 

± 

1.26 

0.74 

± 

0.59 

± 

0.57 


TABLE I. Light quark (uds) Collins asymmetries obtained by fitting the U/L and U/C double ratios as a function of ( 21 ,^ 2 ) 
for kaon pairs in the RF12 frame (upper table) and in the RFO frame (lower table). In the first two columns are reported the 
2 bins and their respective mean values for the kaon in one hemisphere; in the following two columns, the same variables for 
the second kaon are shown; in the fifth column is summarized the ratio of mean values (sin^ dth( 2 ))!(1 + cos^ dth{ 2 )), and the 
asymmetry results are reported in the last two columns. The quoted errors are statistical and systematic, respectively. The 
mean values of the quantities reported in the table are calculated by summing the corresponding values for each KK pair and 
dividing by the number of KK pairs that fall into each ( 21 , 22 ) interval. Note that the and results are strongly 
correlated since they are obtained by using the same data set. 
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[0.15,0.2] 

0.174 

[0.15,0.2] 

0.174 

0.794 

0.19 

± 

0.77 

± 
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± 
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± 
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± 

0.62 

± 
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± 

0.55 

± 

0.23 
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0.174 
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± 
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± 
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± 
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± 
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± 
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± 
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4.12 

± 
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± 

0.35 
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± 

0.63 

± 
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± 
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± 
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± 
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± 
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± 

0.45 

± 
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[0.3, 0.5] 

0.380 

0.789 

4.94 

± 

0.47 

± 

0.39 

2.79 

± 

0.42 

± 

0.24 

[0.2, 0.3] 

0.245 

[0.5, 0.9] 

0.611 

0.782 

7.56 

± 

0.61 

± 

0.52 

4.13 

± 

0.52 

± 

0.32 

[0.3, 0.5] 

0.380 

[0.15,0.2] 

0.174 

0.791 

3.76 

± 

0.55 

± 

0.39 

2.04 

± 

0.50 

± 

0.25 

[0.3, 0.5] 

0.380 

[0.2, 0.3] 
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0.789 

4.24 

± 

0.47 

± 

0.39 

2.44 

± 

0.42 

± 

0.24 
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± 
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± 
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± 

0.45 

± 
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7.70 

± 
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± 
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3.46 

± 

0.57 

± 

0.34 

[0.5, 0.9] 

0.612 

[0.15,0.2] 

0.174 

0.784 

4.75 

± 

0.63 

± 

0.55 

2.39 

± 

0.57 

± 

0.35 

[0.5, 0.9] 

0.611 

[0.2, 0.3] 
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0.782 

8.65 

± 

0.63 

± 

0.52 

4.57 

± 

0.52 

± 

0.32 

[0.5, 0.9] 

0.612 

[0.3, 0.5] 

0.379 

0.781 

9.74 

± 

0.92 

± 

0.57 

5.01 

± 

0.63 

± 

0.34 

[0.5, 0.9] 

0.615 

[0.5, 0.9] 

0.615 

0.777 

12.83 ± 1.54 ± 0.81 

5.06 

± 

0.85 

± 

0.46 

Zl 

(zi) 

22 

(22) 

(sin^ 62 ) 

(10 

-2 

) 

aUC /in 

,-2 


(1 + cos^ 62 ) 





) 

[0.15,0.2] 

0.174 

[0.15,0.2] 

0.174 

0.732 

0.70 

± 

0.78 

± 

0.36 

0.37 

± 

0.69 

± 

0.21 

[0.15,0.2] 

0.174 

[0.2, 0.3] 

0.246 

0.736 

0.48 

± 

0.61 

± 

0.32 

0.23 

± 

0.55 

± 

0.17 

[0.15,0.2] 

0.174 

[0.3, 0.5] 

0.380 

0.748 

1.04 

± 

0.53 

± 

0.31 

0.52 

± 

0.48 

± 

0.17 

[0.15,0.2] 

0.174 

[0.5, 0.9] 

0.612 

0.752 

3.66 

± 

0.64 

± 

0.39 

2.00 

± 

0.56 

± 

0.25 

[0.2, 0.3] 

0.246 

[0.15,0.2] 

0.174 

0.727 

1.68 

± 

0.63 

± 

0.32 

0.93 

± 

0.56 

± 

0.17 

[0.2, 0.3] 

0.245 

[0.2, 0.3] 

0.246 

0.733 

1.47 

± 

0.50 

± 

0.29 

0.79 

± 

0.45 

± 

0.14 

[0.2, 0.3] 

0.245 

[0.3, 0.5] 

0.380 

0.746 

1.89 

± 

0.46 

± 

0.31 

0.98 

± 

0.42 

± 

0.15 

[0.2, 0.3] 

0.245 

[0.5, 0.9] 

0.611 

0.750 

2.47 

± 

0.57 

± 

0.35 

1.26 

± 

0.50 

± 

0.21 

[0.3, 0.5] 

0.380 

[0.15,0.2] 

0.174 

0.724 

1.19 

± 

0.55 

± 

0.31 

0.67 

± 

0.50 

± 

0.17 

[0.3, 0.5] 

0.380 

[0.2, 0.3] 

0.245 

0.731 

1.86 

± 

0.46 

± 

0.31 

1.01 

± 

0.42 

± 

0.15 

[0.3, 0.5] 

0.379 

[0.3, 0.5] 

0.379 

0.744 

2.60 

± 

0.53 

± 

0.30 

1.39 

± 

0.44 

± 

0.15 

[0.3, 0.5] 

0.379 

[0.5, 0.9] 

0.612 

0.749 

2.96 

± 

0.75 

± 

0.37 

1.30 

± 

0.56 

± 

0.22 

[0.5, 0.9] 

0.612 

[0.15,0.2] 

0.174 

0.717 

1.57 

± 

0.63 

± 

0.39 

0.80 

± 

0.57 

± 

0.26 

[0.5, 0.9] 

0.611 

[0.2, 0.3] 

0.245 

0.725 

2.48 

± 

0.57 

± 

0.35 

1.18 

± 

0.50 

± 

0.21 

[0.5, 0.9] 

0.612 

[0.3, 0.5] 

0.379 

0.738 

4.20 

± 

0.86 

± 

0.37 

2.01 

± 

0.60 

± 

0.23 

[0.5, 0.9] 

0.615 

[0.5, 0.9] 

0.615 

0.745 

5.53 

± 

1.43 

± 

0.52 

2.05 

± 

0.81 

± 

0.35 


TABLE II. Light quark [uds) Collins asymmetries obtained by fitting the U/L and U/C double ratios s a function of {zi, Z2) for 
Ktt hadron pairs in the RF12 frame (upper table) and in the RFO frame (lower table). In the first two columns are reported the 2 
bins and their respective mean values for the hadron (K or tt) in one hemisphere; in the following two columns, the same variables 
for the second hadron (K or tt ) are shown; in the fifth column is summarized the ratio of mean values (sin^ ^th(2))/(l+cos^ ^th(2)), 
and the asymmetry results are reported in the last two columns. The quoted errors are statistical and systematic, respectively. 
The mean values of the quantities reported in the table are calculated by summing the corresponding values for each Ktt pair 
and dividing by the number of Ktt pairs that fall into each ( 21 , 22 ) interval. Note that the and results are strongly 
correlated since they are obtained by using the same data set. 
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TTTT sample 


2 l 

(zi) 

22 

(22) 

(sin^ 0 th) 

aUl (-in 

-2 

) 

aVC Mn 

,-2 


(1+cos^ Sth) 


^V2 




112 


) 

[0.15,0.2] 

0.174 

[0.15,0.2] 

0.174 

0.791 

2.64 

± 

0.37 

± 

0.39 

1.25 

± 

0.26 

± 

0.20 

[0.15,0.2] 

0.174 

[0.2, 0.3] 

0.244 

0.789 

3.72 

± 

0.29 

± 

0.40 

1.74 

± 

0.21 

± 

0.19 

[0.15,0.2] 

0.174 

[0.3, 0.5] 

0.378 

0.786 

4.06 

± 

0.24 

± 

0.43 

1.87 

± 

0.19 

± 

0.21 

[0.15,0.2] 

0.174 

[0.5, 0.9] 

0.617 

0.781 

6.26 

± 

0.34 

± 

0.57 

2.80 

± 

0.23 

± 

0.29 

[0.2, 0.3] 

0.244 

[0.15,0.2] 

0.174 

0.789 

3.76 

± 

0.29 

± 

0.40 

1.76 

± 

0.21 

± 

0.19 

[0.2, 0.3] 

0.244 

[0.2, 0.3] 

0.244 

0.788 

4.69 

± 

0.21 

± 

0.41 

2.17 

± 

0.17 

± 

0.20 

[0.2, 0.3] 

0.244 

[0.3, 0.5] 

0.377 

0.785 

4.99 

± 

0.21 

± 

0.44 

2.24 

± 

0.16 

± 

0.21 

[0.2, 0.3] 

0.244 

[0.5, 0.9] 

0.617 

0.780 

8.27 

± 

0.36 

± 

0.58 

3.57 

± 

0.22 

± 

0.28 

[0.3, 0.5] 

0.378 

[0.15,0.2] 

0.174 

0.786 

4.53 

± 

0.25 

± 

0.43 

2.08 

± 

0.19 

± 

0.21 

[0.3, 0.5] 

0.377 

[0.2, 0.3] 

0.244 

0.785 

4.73 

± 

0.21 

± 

0.44 

2.12 

± 

0.16 

± 

0.21 

[0.3, 0.5] 

0.377 

[0.3, 0.5] 

0.377 

0.782 

6.23 

± 

0.33 

± 

0.48 

2.70 

± 

0.19 

± 

0.23 

[0.3, 0.5] 

0.378 

[0.5, 0.9] 

0.619 

0.777 

9.47 

± 

0.59 

± 

0.62 

3.85 

± 

0.29 

± 

0.30 

[0.5, 0.9] 

0.617 

[0.15,0.2] 

0.174 

0.781 

6.58 

± 

0.37 

± 

0.58 

2.94 

± 

0.24 

± 

0.29 

[0.5, 0.9] 

0.617 

[0.2, 0.3] 

0.244 

0.780 

7.45 

± 

0.35 

± 

0.58 

3.21 

± 

0.22 

± 

0.28 

[0.5, 0.9] 

0.619 

[0.3, 0.5] 

0.378 

0.777 

8.77 

± 

0.59 

± 

0.62 

3.55 

± 

0.29 

± 

0.30 

[0.5, 0.9] 

0.622 

[0.5, 0.9] 

0.622 

0.772 

18.46 ± 1.31 ± 0.98 

6.93 

± 

0.54 

± 

0.43 

2 l 

( zi ) 

22 

(22) 

( sin ^ 62 ) 

A ^^ (10 

-2 

) 

aUC /in 

,-2 


(1 + cos ^ 62 ) 





) 

[0.15,0.2] 

0.174 

[0.15,0.2] 

0.174 

0.725 

1.26 

± 

0.31 

± 

0.22 

0.59 

± 

0.24 

± 

0.14 

[0.15,0.2] 

0.174 

[0.2, 0.3] 

0.244 

0.735 

1.66 

± 

0.25 

± 

0.21 

0.78 

± 

0.20 

± 

0.13 

[0.15,0.2] 

0.174 

[0.3, 0.5] 

0.378 

0.744 

1.41 

± 

0.22 

± 

0.22 

0.65 

± 

0.18 

± 

0.13 

[0.15,0.2] 

0.174 

[0.5, 0.9] 

0.617 

0.750 

2.39 

± 

0.30 

± 

0.27 

1.05 

± 

0.22 

± 

0.18 

[0.2, 0.3] 

0.244 

[0.15,0.2] 

0.174 

0.722 

1.52 

± 

0.25 

± 

0.21 

0.71 

± 

0.20 

± 

0.13 

[0.2, 0.3] 

0.244 

[0.2, 0.3] 

0.244 

0.733 

2.12 

± 

0.20 

± 

0.21 

0.98 

± 

0.16 

± 

0.12 

[0.2, 0.3] 

0.244 

[0.3, 0.5] 

0.377 

0.742 

2.13 

± 

0.20 

± 

0.21 

0.96 

± 

0.16 

± 

0.13 

[0.2, 0.3] 

0.244 

[0.5, 0.9] 

0.617 

0.748 

3.03 

± 

0.30 

± 

0.25 

1.31 

± 

0.20 

± 

0.16 

[0.3, 0.5] 

0.378 

[0.15,0.2] 

0.174 

0.718 

1.58 

± 

0.23 

± 

0.22 

0.73 

± 

0.18 

± 

0.13 

[0.3, 0.5] 

0.377 

[0.2, 0.3] 

0.244 

0.729 

2.09 

± 

0.20 

± 

0.21 

0.93 

± 

0.16 

± 

0.13 

[0.3, 0.5] 

0.377 

[0.3, 0.5] 

0.377 

0.738 

2.64 

± 

0.27 

± 

0.22 

1.13 

± 

0.17 

± 

0.13 

[0.3, 0.5] 

0.378 

[0.5, 0.9] 

0.619 

0.746 

3.89 

± 

0.52 

± 

0.28 

1.58 

± 

0.27 

± 

0.18 

[0.5, 0.9] 

0.617 

[0.15,0.2] 

0.174 

0.712 

2.20 

± 

0.30 

± 

0.27 

0.98 

± 

0.22 

± 

0.18 

[0.5, 0.9] 

0.617 

[0.2, 0.3] 

0.244 

0.724 

3.06 

± 

0.31 

± 

0.25 

1.32 

± 

0.21 

± 

0.16 

[0.5, 0.9] 

0.619 

[0.3, 0.5] 

0.378 

0.735 

3.76 

± 

0.51 

± 

0.28 

1.52 

± 

0.27 

± 

0.18 

[0.5, 0.9] 

0.622 

[0.5, 0.9] 

0.622 

0.743 

6.76 

± 

1.06 

± 

0.41 

2.54 

± 

0.46 

± 

0.26 


TABLE III. Light quark (uds) Collins asymmetries obtained by fitting the U/L and U/C double ratios as a function of (21,22) 
for kaon pairs in the RF12 frame (upper table) and in the RFO frame (lower table). In the first two columns are reported 
the 2 bins and their respective mean values for the kaon in one hemisphere; in the following two columns, the same variables 
for the second kaon are shown; in the fifth column is summarized the ratio of mean values (sin^ 64 ^( 2 ))/(l + cos^ dth( 2 )), and 
the asymmetry results are reported in the last two columns. The quoted errors are statistical and systematic, respectively. 
The mean values of the quantities reported in the table are calculated by summing the corresponding values for each tttv pair 
and dividing by the number of tttt pairs that fall into each (21,22) interval. Note that the and results are strongly 
correlated since they are obtained by using the same data set. 
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0 2 4 6 8 10 12 14 16 



0.15 0.2 0.3 0.5 0.2 0.3 0.5 0.20.3 0.5 0.2 0.3 0.5 0.9 

0 . 9 / 0.15 0 . 9 / 0.15 0 . 9 / 0.15 


0.15 0.2 0.3 0.5 0.9 
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0 . 9 / 0.15 0 . 9 / 0.15 0 . 9 / 0.15 


0.15 0.2 0.3 0.5 0.9 


FIG. 5. (color online). Asymmetries measured in the MC sample in RF12 (top) and RFO (bottom) for KK, Kn, and tttt 
pairs. The upper plots show the U/L double ratio, while the lower plots the U/C double ratio. The 16 {zi, Z 2 ) bins are shown 
on the x-axis: in each interval between the dashed lines, zi is chosen in the following ranges: [0.15,0.2], [0.2, 0.3], [0.3, 0.5], 
and [0.5, 0.9], while within each interval the points correspond to the four bins in Z 2 - We subtract these biases from the 
background-corrected asymmetry, and the statistical errors (represented by the bars around the points) are included into the 
systematic uncertainties. 
















































































